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In this third paper the 1902 lie 3 Scale of Temperatures is evaluated both as to its preci- 
sion and its deviations from the thermodynamic Kelvin Scale. Various thermodynamic 
quantities of lie 3 consistent with the 1962 lie 3 Scale are derived and listed. The correction 
to an observed vapor pressure for small amounts of lie 4 is discussed and tabulated. A 
description is given of the method of multiple variable least squares analysis used for deriving 
the final scale equation and for re-analysis of isotherm data. Finally the present status of 
the 1962 He 3 Scale is discussed along with some considerations for the future. 



1. An Evaluation of the 1962 He 3 Scale 



1.1. Fit of the Input 



1961 L.A.S.L. 
Data 



Vapor-Pressure 



One factor in evaluating the 1962 He 3 Scale is the 
fit of the input vapor-pressure data of Part I [46] 
to the Scale. Figure 1 and table 1 of Part II [6] show 
the deviations of the observed data from the final 
scale as T, 2 (P,)-T 58 (P 4 ). The symbol T 62 is the 
temperature on the 1962 He 3 Scale corresponding to 
a He 3 vapor pressure, P 3 , while T 5 $ is the temperature 
on the 1958 He 4 Scale [2] corresponding to the 
experimentally determined isothermal He 4 vapor 
pressure, P 4 . The standard deviation of the data 
from the scale is 0.25 mdeg and the maximum devia- 
tion over the full range is 0.6 mdeg. The data may 
not scatter completely randomly. For example, the 
data points just below 2 °K are all below the 1958 
He 4 Scale and the points just above 2 °K are all 
above the 1958 He 4 Scale. Hence, if one wished to 
obtain the vapor pressure of He 3 which is most 
probably isothermal with a given He 4 vapor pressure, 
a slightly better value (in terms of consistency with 
the observed data) may be obtained by drawing a 
smooth curve through these data points or by using 
direct interpolation equations as discussed in Part I. 
However, the overall fit of the Scale to the input 
data is very satisfactory in comparison with the 
errors of measurement of the two vapor pressures. 



1 Much of the material in this paper III has been included in a review chapter 
[52] on the development of the scale. 

2 Work performed under the auspices of the U.S. Atomic Energy Commission. 



1.2. Fit of the Experimental Thermodynamic Equa- 
tion (ETE) Scale 

Below 0.9 °K the 1962 He 3 Scale was evaluated by 
examining its fit to the ETE Scale described in 
Part II. The ETE Scale is defined by a linear 
equation fitted to an empirical function, F X (P S} T) } 
of He 3 vapor pressure and temperature, 



F x (P,,T) = (a + bT)/R. 



(1) 



As shown by the full line curve in figure 1 of Part II, 
temperatures calculated from eq (1) are in good 
agreement with the 1962 He 3 Scale; nowhere below 
2 °K do the scales differ by more than 0.4 mdeg. 
the 1962 He 3 Scale, denned by eq (9b) of Part II, 
is therefore in effect an ETE scale from 0.2 to 2° and 
an empirical scale above 2 °K. 

The effect of possible errors in the various terms 
of the equation for the ETE scale, over the tempera- 
ture range from 0.2 to 1.0 °K, should also be con- 
sidered in evaluating the 1962 He 3 Scale. The ETE 
scale was obtained from a least squares fit of the 
He 3 -He 4 vapor-pressure intercomparisons between 
0.9 and 2.0 °K to eq (1). The function F X (P 3 ,T) is 
calculated for each experimental P 3 ,T r 58 -observation; 
deviations of this "observed" value of F x (P Zj T) 
from the fitted value, [(a-\-bT 5S )/R], from eq (1) are 
shown as circles in figure 1. The function 
F x (PzjT) is derived in Part II and in ref. [10] as 



F x (Ps,T) = T[(5/2)ln T- 



•ln Pa-Ko 
+UV L )-j x (C sat )+e] (2) 
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Figure 1. — Deviations from the ETE Scale, eq (1). 

The circles are deviations from the fitted equation of the function F x (Pz,Tb%) 
calculated from eq (2) for each input {Pz,Th%) data point. The central cross- 
hatched area is the 95 percent prediction interval [53] for a prediction of F x {Pz,T) 
from eq (1) as calculated from the ± values for a and b given in eqs (3) and (4) . 
The outermost solid curves represent a change in F x {Pz,T) corresponding to a 1 
mdeg change in temperature and are calculated as 0.001 T(d\nPz/dT) sa t. 



where i is the chemical constant including the 
nuclear spin degeneracy; f x (V L ) is an empirical 
power series representing the theoretical and exact 
integral term, 

j(V L ) = (l/RT) £ V L {dP,ldT)^dT 

involving the molar liquid volume, V L ;f x (C^ t ) is an 
equation for the theoretical and exact double integral 

term of eq (3a) of Part II, f(C Mt ) = (l/RT) V AT 
rT> Ji 

I (Csat/T")dT", based on an empirical power 

series fitted between 0.2 and 2.0 °K to data for 
<7 sat , the specific heat of the saturated liquid ; and e is 
the gas imperfection term, eq (4a) of Part II. 

The least squares values for the coefficients a and 
b of eq (1), and their equivalent expressions from 
eqs (6) and (7) of Part II are 



o- ( C sat dT=l7A459± 0.0058 joule mole" 1 (3) 



a=L { 
and 



6 = 5^(1.0°) =9.0098±0.0038 joule mole" 1 deg" 1 . (4) 

The plus-or-rninus values for a and b were found from 
the fit of the data points to the straight line, a+bT. 
Mood [53] has given an equation for the prediction 
interval for a single prediction of F X (P S ,T) from eq 
(1) for any value of T. The ETE Scale is just 
the set of P 3 values obtained by solving eq (2) 
for In P 3 at any given temperature using the pre- 
dicted value of F X (P 3 ,T). Hence a prediction 
interval for F t (P<$,T) can be expressed as an equiva- 



lent prediction interval in P 3 by neglecting the 
contribution of possible errors in the other terms 
on the right hand side of eq (2) . 

The 95 percent prediction interval for a single 
prediction of F X (P S ,T) [with a 95 percent proba- 
bility of containing the statistically "true" value of 
F X (P B T)] is shown as the central cross-hatched 
area in figure 1. The outermost curves in this 
figure show the change in F X (P S , T) corresponding 
to a one millidegree change in temperature. The 
same prediction intervals have been converted to 
equivalent temperature scale errors as listed in 
table 1. 

The effects of random errors of ±3 percent in 
jx (Csat) and e on individual temperatures below 1° 
are also shown in table 1. 

Systematic errors in f x (C sat ), e and in f x (V L ) 
would be compensated between 0.9° and 2.0° by 
the least squares process of fitting the (P 3 ,T 58 ) data 
to eq (1). The scale below 1° would be skewed by 
such errors as is shown in table 1 for assumed 3 
percent systematic errors in all values of C B&t or e . 
Table 1 also shows the systematic errors below 1° 
resulting from use of the approximate empirical 
function, j x (V L ), instead of the graphically integrated 
values of the true thermodynamic function, f(V L ) 
[see table 4 of Part II]. 

Another possible source of systematic error is a 
smooth deviation of the 1958 He 4 Scale from the 
thermodynamic Kelvin Scale. The effect on the 
ETE scale of adding 2 mdeg to each scale input 
temperature is listed in table 1. 

Table 1. Effect of possible systematic or random errors on the 
ETE scale, eq (1) 

Systematic or random changes in the constants a and b of eqs (3) and (4) and of 
temperatures on the ETE scale below 1 °K are listed for the following arbitrary, 
but plausible, cases: 

Case 1. the fit of the input (Pi,T&) data expressed as the 95 percent confidence 
limit [53] for the prediction of a value of In P at any single temperature (see fig. 
i); 

Case 2. random errors of ±3 percent in/* (Cat) ; 

Case 3. random errors of ±3 percent in e, eq (4a) of Part II; 

Case 4. a systematic error of —3 percent in all Csat values, and hence in /z(C 8a t); 

Case 5. a systematic error of -f-3 percent in all values of e; 

Case 6. the difference between the empirical function, j x {Vl), and the numeri- 



cally integrated values of the thermodynamic term, /( Vl) 



-f 



V L (dP/dT) Bat 



dT; and 

Case 7. the increase of all values of input temperatures between 0.9 and 2.0 °K 
by adding 0.002 deg to each T 5 &. 



Case 


Aa 


A6 


AT 




1.0 °K 


0.8 °K 


0.6 °K 


0.4 °K 


0.2 °K 


1 

2 


J 1 mole 
±0.006 


Jjmole deg 
±0. 004 


mdeg 

±0.4 
0.0 

±0.2 
0.2 
0.2 

—0.1 


mdeg 

±0.3 

±0.1 

±0.1 

0.5 

0.4 

—0.1 

1.7 


mdeg 

±0.3 

±0.2 

0.0 

0.8 

0.5 

—0.1 

1.5 


mdeg 
±0.2 
±0.4 
0.0 
1.0 
0.6 
0.0 
1.2 


mdeg 
±0.1 
±0.4 


3 






0.0 


4 

5 

6 


0.069 
0.073 


—0. 060 
—0. 076 


1.0 
0.5 
0.0 


7_ 


0.094 


—0. 026 


0.7 









The overall effect on the 1962 He 3 Scale of random 
errors in the specific heat of liquid He 3 , the virial 
coefficient equations, and the 1958 He 4 Scale might 
amount to as much as three millidegrees. Durieux 
[54] has analyzed the 1962 He 3 Scale for the effects of 
possible errors and has reached similar conclusions. 
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1.3. Fit of the Argonne Laboratory Vapor-Pressure 
Data 

The He 3 -He 4 vapor-pressure intercomparisons of 
Abraham, Osborne, and Weinstock [5] have been 
used in deriving all previous He 3 temperature scales, 
but these were rejected for the 1962 He 3 Scale 
because of apparent thermodynamic inconsistency 
(see refs. [10 and 45]) with recent precise measure- 
ments of the latent heat of vaporization of He 3 by 
the same workers [4, 55]. The fit of the AOW 
vapor-pressure data to the 1962 He 3 Scale is shown 
in figure 2. Above 1.9° the agreement is within 0.5 
mdeg. Below 1.5° the AOW data deviate systemati- 
cally from the 1962 Scale indicating either that their 
He 3 vapor pressures were lower or that their He 4 
pressures were higher than the data of Sydoriak and 
Sherman (Parti). For encircled points in figure 2, 
a He 4 vapor pressure bulb was fastened to the He 3 
bulb. The He 4 pressure may have been high because 
of insufficient correction for effects due to the reflux- 
ing superfluid film. For the uncircled crosses and 
dots, the He 4 bath pressure was measured in the 
pumping line in a warm part of the apparatus. The 
corrected He 4 pressures as published may have been 
high because the effluent gas was significantly colder 
than the tip of the pressure-sensing tube, thus in- 
creasing the actual thermomolecular pressure 
difference. 



1.4. Fit of Heat-of- Vaporization Data 

Heat-of- vaporization, L, data may be used to test 
the thermodynamic consistency of the temperature 
scale. The equations used may be either 



L=T(S G -S L ) 



(5) 



with the vapor entropy, S G , 
Sackur-Tetrode equation, or 



calculated from the 



L=T(V c -V L )(dP/dT)^ 



(6) 



from the Clausius-Clapeyron equation. The use of 
these equations for testing He 4 scales of temperature 
has been discussed extensively; see van Dijk and 
Durieux [25], Durieux [15], Berman and Mate [56], 
and Keller [43]. 

The liquid entropy at 1.0 °K (or any other tem- 
perature) may be calculated from eq (5) and data 

for <7 sat 



as 



S L a.0°)=S o (T)-(L/T)- ( T (C mt /T')dT'. 

J 1-0 



(7) 



The Argonne (WAO) heat-of-vaporization measure- 
ments [4] were undertaken in order to determine 
the entropy of liquid He 3 in this way. 

At one stage of the derivation of the 1962 He 3 
Scale, it was proposed [10] to use all the WAO 
heat-of-vaporization data and eq (7) to get an average 
value for # L (1.0°), which is just the coefficient, 
b, in eqs (1) and (4). Then, values of (a/R) = 



F x (P 3 ,T)—(bT/R) from eq (1) were computed for 
every vapor-pressure datum point between 0.9 and 
2.0 °K. This two-step method failed to converge 
on a stable, consistent pair of constants a and b> 
and led to the decision to undertake the new L.A.S.L. 
intercomparisons of the vapor pressures of He 3 and 
He 4 . Analysis of the new intercomparisons [45] by 
this same two-step method also failed to yield a 
consistent set of a and b values, although these 
data showed much less scatter, as shown in figure 3. 
The heat-of-vaporization data were not used to 
determine the 1962 He 3 Scale; instead a and b 
values were obtained by a least squares fit of the He 3 
vapor-pressure data to eq (1) using for T the T 5S 
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Figure 2. — The deviation in millidegrees of the 1950 Abraham^ 
Osborne, and Weinstock [5] (P 3 , Pa) data front the 1962 He 5 
Scale: T«(P 8 )-T M (P4). 

The uncircled crosses and small dots arc data points below the X-poinl for which 
the Ho 4 cryostat bath pressure was published as Pi. For the encircled points a 
He 4 vapor pressure bulb was fastened to the lie 3 bulb. 



SS DATA 
S L (I.0°) = 1.0720 R 



T,°K 



Figure 3. — Demonstration of the remaining thermodynamic 
inconsistency between eqs (1) and (7) involving the He 3 latent- 
heat data, (L, P3), of Weinstock, Abraham, and Osborne [4] T 
[55]; the (P3, P4) intercomparisons of Sydoriak and Sherman 
[46]; and the 1958 He 4 Scale [2], T 58 (P 4 ). 

The indicated value of Sr,(1.0°) is the average of nine values calculated from eq 
(7) for the WAO (L,Pz) data and is about 1 percent lower than the 1962 He 3 Scale 
value for that entropy, eq (4). The solid circles are the individual values of the 
theoretical constant, (a/R), in eq (1) calculated for each (P3,T5s) data pointlas 
F x (P3,Tm)-[Sl(I.Q°) TmIR] using the indicated average value of SlO-.O ). 
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values that corresponded to the He 4 vapor pressures 
observed concurrently with the He 3 vapor pressures. 
The 1958 WAO heat-of- vaporization data may be 
used to test the 1962 He 3 Scale in the manner ex- 
plained by Durieux [15]. In terms of an apparent 
heat of vaporization, L a (defined as the heat necessary 
to evolve a mole of gas outside of a calorimeter), 
eq (6) becomes 

L a =TV G (dP/dT) SSit 

= -R(l + BIV G +C/V%)[d\nP/d(llT)} bat 

=RT 2 Z(dlnP/dT) sat (8) 

where the compressibility coefficient, Z= (PV / 
RT)=(l+BjV a +C/Vg) from the density virial form 
of the equation of state. 

As shown in figure 4, the observed values of L a 
average about 0.5 percent higher than the values 
calculated from the 1962 He 3 Scale using eq (8), 
a deviation considerably in excess of the estimated 
0.1 percent accuracy [55] of the L a measurements. 
Durieux [15] calculates that the L a values of Berman 
and Mate [56] for He 4 average 0.76 percent higher 
than the values calculated from the 1958 He 4 Scale 
between 2.2 and 3.0 °K. These differences may be 
due in part to a deviation of the 1962 He 3 and 1958 
He 4 Scales from the true thermodynamic temperature 
scale. An estimate of such a difference may be 
obtained by neglecting the variation of (d In P/dT) s&t 
and Z in eq (8). In this case dT/T^ (1/2) (8L a )/L a . 
From the dashed curve drawn through the He 3 data 
of figure 4 this estimate of 8T varies from 3 mdeg 
at 1.2 °K to 1 mdeg at 2.0 °K, a not unreasonable 
range of deviations. For the He 4 data, this approxi- 
mation yields values from 8 mdeg at 2.2 °K up to 
11 mdeg at 3.0 °K which are much larger than ap- 
pears reasonable for the departure of the 1958 He 4 
Scale from the thermodynamic scale. An error 
in the temperature scale also would cause both 
(d In P/dT) B&t and Z to change. 

1.5. Fit of Gas Thermometer, Isotherm, and Acoustic 
Interferometer Measurements 

In principle, gas thermometer, isotherm and acous- 
tic interferometer measurements of the absolute 
temperature associated with a He 4 vapor pressure 
may be used to check the 1962 He 3 Scale by directly 
interpolating the He 4 vapor pressure to an iso- 
thermal He 3 vapor pressure. The direct interpola- 
tion equations and table described in Part I are 
usable for this interpolation of He 4 vapor pressures. 
Since the 1962 He 3 Scale has been made to agree so 
closely with the 1958 He 4 Scale, at least in terms of 
the (P 3 , P 4 ) data of Part I, little really significant new 
information would be expected from this kind of 
comparison. 

Van Dijk [57] has critically reanalyzed all published 
gas thermometer and isotherm measurements and 
compared them with the 1958 He 4 Scale. The 
majority of these measurements between 1.5 and 
3.3 °K seem to indicate that the thermodynamic 
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Figure 4. — Deviations of the values of L a , the apparent latent 
heat of vaporization (see sec. 1.4), observed by Weinstock, 
Abraham, and Osborne [4, 55] from values of L a calculated 
for the 1962 He 3 Scale from eq (8) using eqs (10) and (11). 



temperature may be several millidegrees higher than 
T 58 , although the data scatter over a range of about 
±10 mdeg from T 58 . 

Two preliminary acoustic interferometer results of 
Cataland and Plumb [58] indicate thermodynamic 
temperatures that are 3 ±2 mdeg higher than T 5S 
at 2.0 and 2.2 °K. 

The isotherm measurements of Keller [11, 12] are 
generally conceded to be the most accurate in the 
liquid helium temperature range, and are the only 
data for gaseous He 3 . During the derivation of the 
1962 He 3 Scale, the isotherm data of Keller were 
reanalyzed [10, 48] by Deming's method [47] of least 
squares adjustment with errors in more than one 
measured variable (see appendix A). For the iso- 
therm measurements the observed pressures for each 
data point were corrected or normalized to a cal- 
culated value which would have been observed if the 
cell volume were at the normalized temperature for 
the set of data points. This normalized temperature, 
Tn,58, was taken as the temperature on the 1958 He 4 
Scale corresponding to P 4 , the average He 4 _vapor 
pressure for the set of data. The values of P 4 and 
T N , 58 are listed in tables 2 and 3. The method of 
normalization was essentially the same as that 
described in detail by Keller [11]. 

For the isotherm measurements the quantity 
minimized by the least squares adjustment was 
H[W p (P -P c ) 2 +Wn(N -N c ) 2 } where the independ- 
ently observed variables are P , the normalized 
helium gas pressure, and N , the molar density. 
The fitted function was 



F(Pc, N c ) =P-N C RT{\ +BN e + CN 2 c ) = (9) 

where B and C are the second and third virial coeffi- 
cients, T is the isotherm temperature for a normalized 
set of data, and P c and N c are the calculated or 
adjusted values of P and N. The individual weights, 
W p and W n , were calculated from Keller's assign- 
ment of errors in the various observed quantities. 
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Table 2. Summary of reanalysis of Keller's gaseous He 3 isotherms [12] 

The isotherm designations,jthe data points, and the method of normalization for each isotherm are the same as those used by Keller. Pa is the average He 4 vapor 
pressure for each isotherm and P% is _the average He 3 vapor pressure for the a. isotherm. Tw(P\) is the normalized temperature for each isotherm and is the 1958 He 4 
Scale temperature corresponding to Pi. Pz(P\) is the directly interpolated lie 3 vapor pressure corresponding to Pi, using the direct interpolation equations described 
in Part I, and T 6 2[P3(P4)1 is the 1962 He 3 Scale temperature corresponding i o PziPi). The values of 7 1 i ao , B, and Clisted under the heading "Three-consl ant fit" were 
obtained by solving eq (9) using Deming's method of least squares adjustment with errors in two measured variables (see appendix A). The standard errors and the 
weighted variance are calculated in accordance with the statistical formulas given by Deming [47]. The two constant fit values of 7',.,, and B were obtained by solving 
eq (9) assuming C=0. For the For the y, 8, and e isotherms the variance of the I luce-constant fit is lower and the values of Care greater than the magnitude of their 
statistically assigned errors. Therefore the three constant fit solutions were used for fitting He 3 virial coefficient equations, eqs (10) and (11), and for the temperature 
differences shown in figure 2 of Part II. 



Isotherm 


Units 


7 


8 


a 


6 









No. of points 


18 


18 


11 


10 


8 


P* 

r 58 (P 4 ) 


mm Hg 
mdeg 
mm Hg 
mm Hg 
mdeg 
mdeg 


489. 52 
3779. 1 


178. 09 
2985. 3 
607. 70 


36.16 

2154.2 

197. 43 

197.62 

2154. 4±0. 2 

2.1550 


13.070 
1813.4 
105.33 


3.782 
1510.3 


Pz(Pi) 


52.21 


Pi 






T 62 [Ps(Pi)] 




2985. 2±0. 2 


1813. 1±0. 2 


1510. 2±0. 2 


T 62 (P3) 

















THREE-CONSTANT FIT 



Tie, 

B. 
C. 



Variance 

(T 5 8— Ti so ) 

(T 62 -T iB o) 

(T 62 (P3)-Ti S o)-. 



mdeg 

cm 3 /mole 

cm 8 /mole 2 



mdeg 
mdeg 
mdeg 



3782. 3 ±2. 2 
-65. 50±0. 86 
11 92 ±306 
7.0 
-3.3±2.2 



2992. 4 ±1.6 
-86. 31 ±0.52 
1746 ±136 
13.8 

-7.1±1.6 
-7.2±1.8 



2153.7±1.0 
-118.81±1.35 
548±822 
6.1 

0.5±1.0 
0.7±1.2 
1.3±1.0 



1816. 4±0. 8 

-145. 08± 1.61 

2611±1289 

0.86 

-3. 0±0. 8 

-3.2±1.0 



1510. 4±2. 5 

-162.15±10.7 

-11, 069 ±16, 390 

4.8 

-0.1 ±2. 5 

-0.2±2.7 



TWO-CONSTANT FIT 



Tiso 

B 

Variance. 



mdeg 
cm 3 /mole 



3774.4±1.3 
-62.18±0. 19 
13.2 



2974. 4 ±2. 8 
-79.66±0.33 
156 



2153. 2±0. 6 
-11 7. 92 ±0.26 

5.7 



1814.9±0.4 
-141. S7 ±0.29 
1.20 



1512. 0±0. 7 
-169.30d=1.15 

4.4 



Table 3. Summary of reanalysis of Keller's gaseous He 4 isotherms [11] 

The various quantities have the same meaning as in table 2. The three constant fits gave a statistically better fit for only one isotherm, so the results of the two 
constant fits have been accepted as being more meaningful and have been used for figure 2 of Part II. 



Isotherm 


Units 


A 


B 


C 


D 


E 


No. of points 


12 


14 


12 


17 


9 


P 4 


mm Jig 

mdeg 
mm Ilg 
mdeg 


588. 55 
3954. 2 


291.31 
3339. 2 


147. 55 
2864. 

528. 42 
2864. 0±0. 2 


52. 758 
2315.1 
254. 50 
2315. 0±0. 2 


35. 541 


T 5 8(Pl)- 


2147.3 


Pi(pi) 


195.21 


T 62 [P3(P4)] 






2147. 6±0. 2 











TWO-CONSTANT FIT 



T i80 


mdeg 
cm 3 /mole 


3953. 7±2. 1 
-83.97±0.00 
8.7 
-1.5±2.1 


3338.4±2.1 
-103.60±0.78 
17.8 
0.8±2.1 


2863. 5±0. 9 
-124.16±0.63 
3.6 

0. 5±0. 9 
0.5±1.1 


2316. 7±2. 3 

-153.38±4.21 
23.4 
-1.6±2.3 

-1.7±2.5 


2148. 9±0. 9 


B 


-176. 04 ±2. 04 


Variance _. __ __ ______ 


1.69 


(Tss-Tiso) 


mdeg 
mdeg 


-1.6±0.9 


(T 6 2-Ti 80 ) 


-1.3±1.1 











THREE-CONSTANT FIT 



7 iso 

B 

C_ 

Variance _ 



mdeg 

cni :1 /molc 

cm 6 /mole 2 



3953. 8±4. 3 
-82. 16±3.52 
-1452±2773 
9.4 



3335. 3±3. 7 
-99.80±3.78 
-31 33 ±3042 

17.7 



2860. 9 ±1.5 
-118.47±2.86 
-6741 ±3331 

2.8 



2311.9±6.1 
-132. 1 ±25. 7 
50, 490_fc60, 020 
23.9 



2146.5d=2.8 

-162.6±15.8 

-38,900±45,160 
1.75 
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The results of the reanalysis of Keller's gaseous 
He 3 data [12] were fitted to empirical equations for 
the He 3 second and third virial coefficients. These 
equations are 



and 



B= (4.942 -270.976/T) cm 3 /mole (10) 

C= (2866/VT) cm 6 /mole 2 . (11) 



The deviations of both He 3 and He 4 isotherm 
temperatures from the 1962 He 3 and 1958 He 4 
Scales are shown in figure 2 of Part II. Keller 
used a He 4 vapor-pressure thermometer for all but 
one of his isotherms, but his observed P 4 's have 
been interpolated to ZYs by use of the direct in- 
terpolation equations of Part I. The weighted 
average of the isotherm temperatures is 1.50 ±1.0 
mdeg above the corresponding T 58 values and 1.52 ± 
1.2 mdeg above the jT 62 values. For one isotherm 
He 3 containing about 0.25 percent He 4 was con- 
densed in the vapor-pressure bulb. The corrected 
average pressure for pure He 3 was 197.62 ±0.02 mm 
Hg (T 62 =2.1550 °K) at an isotherm temperature of 
2.1537±0.0010 °K. It was concluded that nothing 
would be gained by trying to base a He 3 scale 
more directly on these isotherm temperatures since 
the 1958 He 4 Scale is based on these data. More- 
over, the 1962 He 3 Scale, as it has been set up to 
agree with the 1958 He 4 Scale, adequately expresses 
the experimental He 3 — He 4 vapor-pressure relation 
for most practical purposes. 



Thermodynamic Properties of He 3 
sistent with the 1962 He 3 Scale 



Con- 



Two temperatures which are frequently given in 
discussing helium vapor pressure scales are the 
boiling point and the critical point. The boiling 
point of He 3 on the 1962 He 3 Scale, at a vapor 
pressure of 760 mm Hg at °C and standard gravity, 
is 3.1905 °K. The temperature at the critical 
point for He 3 was measured by Sydoriak and Sher- 
man [46] by observing the He 4 vapor pressure. 



This measured critical temperature on the 1958 
He 4 Scale was 3.3240 ±0.0018 °K at a He 3 critical 
pressure, P c , of 873. 0± 1.5 mm Hg. The calculated 
critical temperature on the 1962 He 3 Scale is the 
temperature given by eq (9b) of Part II [6] corre- 
sponding to P c and is 3.3246 ±0.0017 °K. This 
close agreement with the directly measured tem- 
perature is quite gratifying and somewhat unexpected 
since the scale-defining equation was fitted over the 
full range of the vapor-pressure intercomparisons. 

Table 4 gives a number of quantities which are 
consistent with the 1962 He 3 Scale, including, for 
comparison, values of He 4 vapor pressures from the 
1958 He 4 Scale. The concentration derivative of 
In Px, (d In P x /dX) Tt x=i, where X is the mole 
fraction of He 3 and P x is the corresponding mixture 
vapor pressure, is useful for correcting for the 
He 4 impurity in a He 3 vapor-pressure thermometer. 
This correction is discussed in section 3. 

The table values for (7 sat , the specific heat of the 
saturated liquid, were calculated from the empirical 
equation, eq (3b) of Part II, fitted to specific-heat 
data for the derivation of the 1962 He 3 Scale. 
The calculated values agree with the experimental 
data to ±1 percent; the deviations are given in 
table 2 of Part II. Values for S L , the entropy 
of the saturated liquid, were calculated from the 
relation 

S L (T)=S L (1.0°)+ f (C s JT')dT (12) 

Jl.O 

using the above-mentioned C\ &% equation and the 
1962 He 3 Scale value of S L (1.0°), given in eq (4). 

Values for the second and third virial coefficient 
were calculated from eqs (10) and (11) while V G , 
the molar volume of the saturated vapor, was 
calculated from the equation given for Z in eq (8). 
These calculated values of V G agree to within 
1 percent with the smoothed fit to the experimental 
data of Kerr [59] up to 2.8 °K. The values of 
V L , the molar volume of the saturated liquid, are 
taken from the table of Kerr and Taylor [60]. 
The values of L, the heat of vaporization, are calcu- 
lated from eq (6). 



Table 4. Thermodynamic properties of He 3 and vapor pressures of He 4 consistent with the 1962 He 3 Scale 



T 


Ps 


P4(T 58 ) 


/dluP2,\ 
\ dT ) 


fd In P x \ 

K dX J T 

at 

X=l 

(see notef) 


C Ba t 


Sl 


Virial coefl. 


Vg 


Vl 


L 




B 


C 




°K 




microns Hg at 


°C andstd.g 


deg- 1 




0.23 

0.47 

0.70 

0.88 

0.92 

0.92 

0.91 

0.89 

0.88 

0.86 

0.810 

0.781 

0.754 

0.728 

0.705 


J/deg mole 


J/deg mole 


3.703 
5.731 

7. 071 
8.117 
9.010 
9.824 
10. 605 
11.378 
12.159 
12.951 


cm^/mole 


cmPlmole 2 


cm 3 /mole 


cm s /mole 
36. 8346 
36. 7809 
36. 7203 
36.7197 
36. 7849 
36. 9C73 
37. 0897 
37.3517 
37.7151 
38.201 
38. 854 
39. 680 
40. 734 
42. 124 
44. 049 
46.818 


J/mole 
20. 56 


0.2 


0. 0121 
28.115 
544. 49 
2892. 5 
8842. 4 
20163 
38516 
65467 
102516 
151112 
212673 
288613 
380383 
489549 
617907 




73.32 
21.04 
10.57 
6.650 
4.721 
3.613 
2.908 
2.425 
2.077 
1.815 
1.611 
1.448 
1.317 
1.210 
1.122 


2.728 
3.146 
3.477 
3.812 
4.222 
4.753 
5.414 
6.201 
7.079 
7.983 


-1350 

-672. 5 

-446. 7 

-233. 8 

-266. 

-220.9 

-188.6 

-164.4 

-145.6 

-130.6 

-118.2 

-108.0 

-99.3 

-91.8 

-85.4 


6409 
4532 
37G0 
3204 
2866 
2616 
2422 
2266 
2136 
2027 
1932 
1850 
1777 
1713 
1655 


103X107 

887X103 
68274 
16909 
6776. 6 
3476. 7 
2060. 7 
1338.9 

925.5 

667.8 

496.6 

376.7 

288.8 

221.0 

163.1 


24. 39 


0.4 




27.97 


0.6 
0.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 
2.2 


0. 2812 
11.445 
120.00 
625. 02 
2155.4 
5689. 9 
12466 
23767 
40466 
63304 
93733 
132952 
182073 


31.42 
34.61 
37.51 
40.08 
42. 29 
44.07 
45.34 
45.99 
45.91 


2.4 






44.94 


2.6 






42.84 


2.8 






39.11 


3.0 






32.25 











fHere X is the He 3 mole fraction. The He 4 mole fraction, (1— X), is used to make the impurity correction. 
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The heat of vaporization at absolute zero, L , 
can be obtained from the fitted constant, a, of eq 
(3) if the specific heat integral is evaluated. The 
specific heat of liquid He 3 under a pressure of a 
few centimeters of mercury has been measured by 
Anderson, Wheat-ley, and coworkers [01 J down to 
a temperature of 0.015 °K. Their empirical equation 
for C p , the specific heat at a constant pressure of 
0.12 atm, is (C P /E)=2.89T-7.80T 2 +7.09T 3 for 
0<r<0.3°. Integration of this equation yields 

r*C m 4T** f ' 2 CV/r=0.331=L0.013 J/mole. This 
Jo # ^ Jo 

value is in excellent agreement with the value 
0.333 ±0.010 J/mole obtained by graphical inte- 
gration of the smooth curve through recent specific 
heat data given by Strongin, Zimmerman, and 
Fairbank [62]. From the empirical equation for 

/•l.O 

C B&t used for the T 62 scale, (7 Bat d2 n, = 2.784±0.056 

JO. 2 

«7/mole. Our corresponding value of L is 20.56 
±0.07 J/mole. 
The entropy of ( he saturated liquid can be assumed 

to be approximately equal to {C v \T f )dT f from 

Jo 
the 0.12 atm C v equation of Anderson et al. [61] 
given above. From this approximation 5^(0.23 °K) 
is 4.05±0.17 J/mole (leg. Strongin et al. [62] 
also have computed the entropy from their extrap- 
olation of specific heat data to °K. Their value 
for the entropy at 0.23 °K is 3.97±0.17 J/mole deg. 
These values are in good agreement with the value 
4.02 ±0.25 J/mole deg obtained by Weinstock, 
Abraham and Osborne [4], who used eq (5) to cal- 
culate S L (\.5°) from their three separate latent 
heat measurements at that temperature. The 
equivalent of eq (12) was used by WAO to compute 
other entropies in the range of the then available 
specific heal data. The value of #7,(0.23°) con- 
sistent with the 1962 He 3 Scale, calculated from eq 
(12),isS L (Q.23°) = 4.Q9±().]() J/mole (leg, estimating 
a ±2 percent limit of error for the specific heat 
integral term. The good agreement among these 
values is gratifying but it may be fortuitous. For 
example, the He 3 temperature scale used by 
Weinstock et al. [4 J to assign temperatures to 
observed He 3 vapor pressures was the T L Scale, 
eq (9) of ref. [9]. The deviations of the T L Scale 
from the 1962 He 3 Scale are shown in figure 5 and 
range from ±6 to —3 mdeg. At 1.5 °Iv, however, 
the two scales are in good agreement as to the 
temperature. Thus the value of J 7 used in computing 
both S G and L/T in eq (7) would give essentially 
the same value for liquid entropy at T^1.5 °K on 
the T L and 1962 He 3 Scales. If, however, a value 
of L calculated from eq (6) is used in eq (7) instead 
of a WAO experimental value of L, significantly 
different values are obtained for S L depending 
upon which of the two scales, T L or 1962 He 3 , is 
used for the calculation of (dPs/dT) B&t since there 
is a difference of about 0.4 percent in these derivatives 
on the two scales as inferred from the slope of the 
AT curve in figure 5 (see, for example, ref. [63]). 
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Figure 5. — Graphs of differences between old and new He 3 and 
He 4 temperature scales. 

The Be 3 curve, TL{Pi) — T^(P^), is a graph of the difference between the Tl 
Scale [9] and the 1962 Be 3 Scale. The lie 1 curve is I he difference between the 
Tlis Scale [25] and the 1958 He 4 Scale given in table 6 of ref [2]. It was intended 
in their derivations that t lie '/'/. Scale would reproduce t he 7 '/.v. Scale, and that the 
1962 Bea Scale would reproduce the L958 Be* Scale. The difference between the 
two curves is explained principally by the new (7'./':) intercoinparison data of 
Sydoriak and Sherman [46] used in deriving the 1962 Be 3 Scale. 

3. Corrections to the Measured Pressure of 
a He' Vapor-Pressure Thermometer for 
the He 1 Impurity 

Sydoriak and Sherman have discussed in Part I 
corrections to be applied to an observed pressure 
as measured in the laboratory in order to obtain 
the vapor pressure at the surface 4 of liquid He 3 . 
The correction for a small He 4 impurity will be 
discussed in more detail here. 

The presence of He 4 in liquid He 3 thermometers 
lowers the vapor pressure below that of pure He 3 . 
Much of the He 3 available for purchase up to the 
present has contained significant quantities of He 4 ; 
even the He 3 used for the 1962 He 3 Scale input 
data measurements [46] contained about 0.04 per- 
cent He 4 . Specially purified He 3 containing less 
than 0.01 percent is being made available for 
purchase through the Division of Research of the 
United States Atomic Energy Commission [64]. 

The correction for small amounts of He 4 may be 
calculated from the approximate relation 



>(l-X)(dP x /dX) 



(13) 



where P x is the observed vapor pressure for mole 
fraction, X, of He 3 . Although very little data 
has been published for the vapor pressures of dilute 
solutions of He 4 in He 3 , the liquid phase diagrams 
of Sydoriak and Roberts [32] and EsePson and 
Berezniak [38] indicate that (dP x /dX) T is probably 
a constant for X>0.9 and hence eq (13) may be 
valid for concentrations of up to 10 percent He 4 . 
Smooth values of this derivative in reduced form, 



(d In P x /dX) TtX=1 =(l/P 3 )(dP x /dX) T , x = l 

~[Hr)y- X) > (I4) 



573 



are listed in table 4 and are compared with the raw 
data of Sy doriak and Roberts [32] in figure 6. 
Between 0.6 and 2.0 °K the table 4 entries and the 
solid curve in figure 6 were calculated from eq (14) 
using the X=0.9 values from the smoothed table 
of ref. [32]. Below 0.6 °K the derivative has been 
assumed to go linearly to zero at °K mainly because 
such an extrapolation is consistent with the data 
below 1 °K, as shown in figure 6. Actually for all 
practical purposes this derivative may be zero 
below about 0.2 °K since at these temperatures 
liquid He 3 — He 4 mixtures with very small amounts 
of He 4 undergo phase separation [65] and, when 
this happens, the vapor pressure becomes independ- 
ent of further decreases in X. Above 2.0 °K, 
the table 4 values were calculated from Raoult's law, 



(dlnPJdX) T , x=l =l-(P 4 /Pz), 



(15) 



shown as the dashed curve in figure 6. That law 
is in fairly good agreement with all existing vapor- 
pressure data for X>0.89 and T> 1.7 °K. 

The concentration of He 4 in the saturated vapor 
will be much less than that of the liquid. This 
fact can cause fractionation and may require further 
correction in precise thermometry if much of the 
He 3 in the thermometer system is not condensed. 

The correction to a temperature measurement 
made with a He 3 sample containing 0.1 percent 
He 4 in the liquid phase ranges from 0.02 mdeg at 
0.4 °K to 0.71 mdeg at 3.2 °K. For most purposes 
the correction may be calculated from the approxi- 
mate relation 

P 3 -P*«(l-Z) P x (d In P x /dX) TtX=l . (16) 

The derivative below 2 °K may be interpolated 
from a curve through the values listed in table 4 
or from figure 6. Above 2 °K the derivative can 
be easily calculated from eq (15). A careful experi- 
mental determination of the correction above 2 °K 
would improve the accuracy of this correction. 




Figure 6. — The correction factor for He 4 impurity. 

The circles are the data points of Sydoriak and Roberts [32] for He 3 mole frac- 
tion, X, of 0.897. The solid curve from 0.6 to 2.0 °K is drawn through the 
smoothed table values of ref. [32] for X=0.9. The long and short dashed curve 
is an arbitrary extrapolation to ° K. The dashed curve is Raoult's Law, eq (15) . 
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Figure 7. — Correction for He 4 impurity, (T 6 2 — T u ), to be added 
to the uncorrected temperatures, T u , calculated from the ob- 
served vapor pressure and the 1962 He 3 Scale. 
The correction is calculated as T 6 2-Tu=(l-X)(dlnP x /dA')r,.Y-i/(dlnP 3 /dT) 

based on eqs (13) and (14) and using the derivative values listed in table 4. The 

curves can be used for ten times the indicated percent of He 4 in the liquid if the 

ordinate is multiplied by ten. 



Figure 7 may be used to read the temperature 
correction directly for known He 4 concentrations 
in the liquid. 



Status of the 1962 He 3 
peratures 



Scale of Tem- 



The 1962 He 3 Scale was presented [66] at the 
Eighth International Conference on Low T Tem- 
perature Physics in London, September 17-21, 1962. 
During the next week the new scale was discussed in 
detail by the Advisory Committee on Thermometry 
of the International Committee on Weights and 
Measures. A scale proposal along with ref. [54] 
has been published along with the minutes of the 
meeting of the Advisory Committee in Sevres, 
France, September 26 and 27, 1962. Quoting from 
the draft of those minutes [67]: 

"II a estime que V Echelle 3 He 1962 doit egalement 
etre recommandee pour Tusage general, avec la 
designation T 62 . 

"Les deux echelles T 58 (VEchelle 'He 1958) et T 62 
peuvent etre utilisees concurremment dansle domaine 
ou elles sont valables. Cependant, quand il s'agit 
de Tadoption de cette nouvelle echelle 3 He comme 
partie de TE.I.P.T., on doit prendre soin d'eviter 
toute ambiguite dans le domaine de recouvrement 
avec F echelle 4 He." 3 

The recommendation by the Advisory Committee 
of the 1962 He 3 Scale was approved in October 1962, 
by the International Committee on Weights and 



s "It (the Advisory Committee) has deemed that the 1962 He 3 Scale ought 
equally to be recommended for general usage with designation Tw. 

"The two scales Tm and T&2 can be used concurrently in the range where they 
are valid. However, when it is a question of adopting this new He 3 scale as a 
part of the I.P.T.S. (International Practical Temperature Scale), care should be 
taken to avoid any ambiguity in the region of overlap with the He 4 scale. 
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Measures meeting in Sevres. The International 
Committee requested the United States and Russian 
Governments to take steps to make high-purity He 3 
available internationally for vapor-pressure thermom- 
etry, and to prepare and distribute known mixtures 
of He 3 and He 4 for the calibration of or the checking 
of the calibration of apparatus for measuring the 
isotopic purity of He 3 . t 

The status of the 1962 He 3 Scale and proposed 
changes in the International Practical Temperature 
Scale have been discussed by Stimson [68] and 
Brickwedde [69]. 

5. Practical Lower Limits for He 3 and He 4 
Vapor-Pressure Thermometry 

Practical lower limits for helium vapor-pressure 
thermometry are determined by the desired precision 
of temperature determination, the accuracy with 
which the vapor pressure can be measured, and by 
how small or how well known the temperature 
difference is between the thermometric fluid and the 
object of the measurement. 

Below about 2.2 °K, He 3 can much more easily be 
made isothermal with an object the temperature of 
which is being measured than can He 4 , for the reasons 
discussed in Part I. A major factor in fixing- a 
practical lower limit for He 3 thermometry is the time 
necessary to achieve an equilibrium thermomolecular 
ratio, PJPw, of cold pressure, P c , to warm pressure, 
P w , and more importantly, of sufficient data to 
permit accurate calculation of that ratio. 

To measure temperatures from 0.3 to 0.25 °K to a 
millidegree precision using He 3 , vapor pressures from 
1.9 to 0.24 jj, would have to be measured to 4 or 5 
percent accuracy. This accuracy in the pressure 
system is relatively easy to achieve. However, the 
error in PJPw may amount to several percent because 
of marked quantum effects in the gas viscosity and 
because no He 3 measurements have been made of the 
empirical thermal transpiration coefficients postu- 
lated by Weber [70]. The He 3 thermomolecular 
pressure measurements made by Roberts and 
Sydoriak [37] showed that for a cold temperature of 
2 °K observed He 3 P c /P„ ratios from 1.0 to 0.5 
agreed within 1 percent with ratios calculated for 
He 4 using the Weber-Schmidt equation [71]. 

Roberts and Sydoriak also discuss the possible 
effect of the low-temperature viscosity on pressure 
ratios in the tubing at liquid helium temperatures. 
For warm and cold temperatures, T w and T c , below 
5 °K and pressure ratios close to 1, eq (5) of ref. [37] 
is (pjp w ) 2 = 1-0.005642 (TJ-T c 2 )/(Rp w ) 2 where 
R is the tubing radius in centimeters and p w is in 
microns of mercury. In order to be reasonably sure 
that the calculated thermomolecular ratio is within 
4 or 5 percent of the true ratio for He 3 it is suggested 
that the calculated (p e /Pv>) be kept greater than 0.5 
in the rest of the pressure sensing tube. For the 
stepped pressure sensing tube used by Sydoriak and 
Sherman with dimensions shown in figure 4 of Part I 
[46], this suggests a practical lower limit of Rp c = 
0.155 cm microns or a vapor pressure of 1.16 \x 
(T 62 = 0.287°)., 



The practical lower limit for He 4 vapor-pressure 
thermometry is strongly dependent on numerous 
details of the particular technique and apparatus 
used. In many cases experimental temperatures are 
determined from the vapor pressure of the bath in 
which the experimental apparatus is immersed. For 
bath temperatures and with proper corrections for 
pressure drops associated with film reflux and/or 
motion of effluent vapor from the bath, temperatures 
may be determined down to about 1 °K. Figure 1 
of Part I shows errors due to film reflux pressure 
drops for typical sizes of pressure sensing tubing. 
However, for temperatures determined by thermal 
contact with a closed-bulb He- vapor-pressure ther- 
mometer the practical lower limit is raised to just 
below the A point, because of the large and variable 
temperature discontinuity at the bulb wall associ- 
ated with recondensatiorj of the refluxing film. 
Figure 1 of Part I also shows this "correction" for 
a variety of conditions., 

If the bulb wall is a poor thermal conductor, a 
still higher limit is imposed due to the effects of the 
density maximum at about 2.18 °K, discussed by 
Chase, Maxwell, and coworkers [72]. 

In the case of precision measurements using He 3 
thermometry, the publication of only a few pertinent 
details, such as the diameters and temperatures of 
any stepped junctions in the pressure sensing tube, 
will enable future correction of data as improved 
thermal transpiration corrections become possible. 

6. Conclusion and Considerations for the 
Future 

These papers on the development of the 1962 He 3 
Temperature Scale have emphasized marked differ- 
ences between liquid He 3 and He 4 . One is in the 
low-temperature variation of the liquid entropy. 
The quantum mechanical exchange forces in liquid 
He 3 lead to appreciable correlation of the nuclear 
spins in the liquid at temperatures as high as 1 °K. 
The difficulty in fitting the entropy curve with a 
simple analytical expression has led to difficulty in 
attaching theoretical significance to certain terms in 
empirical vapor-pressure equations at temperatures 
below the temperature of fitted data. 

The superfluid transition of liquid He 4 makes 
accurate isothermal intercomparisons of He 3 and He 4 
vapor pressures difficult below the He 4 X-point. The 
1961 intercomparison data of Sydoriak and Sherman 
[46], however, are considered to be valid to a few 
tenths of a millidegree. The 1962 He 3 Scale, which 
is based on these intercomparisons, is believed to be 
in agreement with the 1958 He 4 Scale to within two 
or three tenths of a millidegree over the full range of 
the intercomparison data, from 0.9 °K to the critical 
point, 3.324 °K. A lower limit for practical He 3 
vapor-pressure thermometry is 0.25 °K corresponding 
to a He 3 vapor pressure of 0.24 n Hg. 

A number of independent measurements, including 
measurements of the latent heats of vaporization of 
He 3 and He 4 , precise isotherm measurements, and 
preliminaiy acoustical interferometer measurements, 
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indicate that both the 1962 He 3 Scale and the 1958 
He 4 Scale may be 2 or 3 mdeg lower than true 
thermodynamic temperatures between 1 and 3.3 °K. 
When more accurate data on the thermodynamic 
properties of He 3 and He 4 are available, together 
with better determinations of absolute thermo- 
dynamic temperatures, the precise intercomparison 
data will allow both the He 3 and the He 4 scales to be 
adjusted consistently and simultaneously. He 3 and 
He 4 scales, thermodynamically consistent between 
1 and 3 °K, may be extended by precision para- 
magnetic salt thermometry above 3 °K using He 4 
and below 1 °K using He 3 . Also, the scales may be 
extended downward by thermodynamic calculations, 
using coefficients evaluated between 1 and 3 °K as 
in the development of the 1962 He 3 Scale and the 
T L55 He 4 Scale [25]. Further measurements on He 3 
that would be helpful in improving the Scale are 
measurements of the vapor pressure and thermo- 
molecular pressure ratio down to temperatures of 
0.25 °K, more accurate specific heat measurements 
on the liquid below 2 °K, and specific heat, latent 
heat, and high-precision vapor-density measurements 
above 2 °K. 

7. Appendix A. Method of Multiple Vari- 
able Least Squares Analysis 

In the usual least squares adjustment of experi- 
mental data the implicit assumption is made that 
only one measured quantity, often called the de- 
pendent variable, is in error. This dependent vari- 
able is measured as a function of one or more inde- 
pendent variables. The usual procedures for the 
adjusting of sets of observations to an empirical or 
theoretical equation assume no error in the inde- 
pendent variables. Frequently, however, the "in- 
dependent variable " is also determined by a physical 
measurement and is just as subject to error as the 
dependent variable. 

Deming [47] has discussed a method of adjust- 
ment of sets of observations of two or more variables 
when all variables are subject to error. His method 
is essentially equivalent to statistical weighting of 
the data. The method has been coded by P. 
McWilliams of the Los Alamos laboratory for use 
with the Gaussian method of iterative least squares 
analysis adapted by Moore and Zeigler [73] for high- 
speed computers. We will give a simple example 
of the method of analysis with exaggerated scatter 
and uncertainty. 

7.1. Simple Example 

As a simple illustration consider fitting the best 
straight line to the following three points: 



X 


Y 


1±1 


1±1 


4±1 


5±1 


5il 


4±1 




Figure 8. — Simple example of the least squares Jit of a straight 
line to the three observed points shown as dots. 

The dashed line is the usual least squares solution of Y=a-\-bX which mini- 
mizes S(Fo-y c ) 2 . The solid line is the solution assuming both X and Fare 
subject to equal absolute errors. 



The dashed line in figure 8 shows the fit of the usual 
least squares solution of Y=a-\-bX assuming that 
only Y is in error and equally weighting the points. 
The solid line, Y=X, is the solution assuming that 
the X and Y values are all equally subject to equal 
absolute errors. The general case of fitting to a 
straight line is discussed in detail by Deming. For 
our special case of equal weighting, the two solutions 
go through the centroid but have different slopes. 

Actually the quantity which is minimized by the 
least squares solution is 

M 
S = J2 [W*i<Xo,i-Y e ,d'+W Xti (Xa. i - 7 g< ti y] (17) 



where the subscript o indicates the observed value; 
subscript c indicates the calculated or adjusted value; 
and w Vt i and w Xt t are the weights of the ith of M sets 
of observations. Thus in figure 8 the solid line is 
identified as minimizing 

S=S[(AF) 2 +(AX) 2 ] 

since all the individual weights are equal. AY is 
defined as (Y -Y c ) and AX=(X -X C ). 

The geometrical relation between AY and AX 
that is assumed in Deming's method of solution is 
illustrated in figure 18 of ref. [47]. The relation for 
each point depends on the individual weights for 
that point. 

If w x =w y , this condition requires that the line 
from the observed point, (X,F ), to the calculated 
point, (X CJ Y C ) } be perpendicular to the calculated 
curve. 



576 



7.2. General Discussion of Method 

In many cases the least squares problem can be 
expressed as minimizing eq (17) subject to M equa- 
tions: 

F(X eii ,Y Cii )=0. (18) 

As an example consider a case in which three coeffi- 
cients, a, b, and c, are to be fitted by least squares to 
the relation Y=f(X, a, b, c). We take 



F(X,Y)=Y-j(X,a,b,c). 



(19) 



Using his geometrical relation between AY and AX, 
and using the Gaussian method of linearization of 
the problem by a truncated Taylor's Series expansion 
of the function, Deming derives a typical least squares 
solution "normal" equation as 

m m 

A<£ W t Fl i+Abj: W t F a , % F C , t 

M M 

+Ac J2 W t F a , t F c , 4 =S W t F a , ( F , ,. (20) 

i=l i=l 

In this equation, the subscripts a, b, and c designate 
the partial derivatives of F with respect to the 
respective coefficients; Aa=ai—a c , the difference 
between the initial value, a u and the calculated or 
adjusted value, a c ; F 0ti =F(X 0ti , Y 0ti , a u b u c 2 ); 
and the weights are 



W i =(F 2 x , t lw x , t +Fl. t /w y . t y 



(21) 



F Xti and F Vfi are the partial derivatives of F with 
respect to X and F, respectively. Where a set of 
observations, X oA and Y 0ti , can be assigned re- 
spective standard deviations a Xti and a yJ , we take 
the weights as: 



W Xti =l/<r 2 Xt i] Wy t i=l/cr 2 yi i. 



(22) 



If variations in X and Y are independent, substitu- 
tion of eq (22) in eq (21) yields 



W ( =llrt, t 



(23) 



and hence the method is equivalent to statistical 
weighting of the individual sets of observations. 

Using the Gaussian method, several iterations 
using the successively improved values of the coeffi- 
cients may be used. The partial derivatives are 
usually evaluated using the observed values of X 
and Y in each iteration. 

As pointed out by Deming, an advantage of the 
method is that the solution is unequivocal; i.e., the 
coefficients are independent of variations in the 
way of writing the equation. Thus a straight line 
can be written as either Y=a J rbX or X=— a/b-{- 
Y/b and correspondingly F would be (Y—a—bX) or 
(X-\-a/b—Y/b). The values of a and b and the 
adjusted points will be the same to within higher 
powers of the residuals, AX and AY. 



The least squares code used computes standard 
deviations of the fitted parameters based on the 
fit of the observed points to the calculated points. 
The code also permits computation of standard 
deviations of points calculated at interpolated or 
extrapolated points as well as at data points. 



We are indebted to K. K. Zeigler and P. McWil- 
liams for many discussions and for the development 
of the computer code. 
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